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We encode the many-body wavefunction of a Bose-Einstein condensate (BEC) in the N- 
particle sector of an extended catalytic state. This catalytic state is a coherent state for 
the condensate mode and an arbitrary state for the modes orthogonal to the condensate 
mode. Going to a time-dependent interaction picture where the state of the condensate 
mode is displaced to the vacuum, we can organize the effective Hamiltonian by powers of 
Requiring the terms of order to vanish gives the Gross-Pitaevskii equation. 

Going to the next order, N^, we derive equations for the number-conserving Bogoliubov 
approximation, first given by Castin and Dum [Phys. Rev. A 57, 3008 (1998)]. In contrast 
to other approaches, ours is well suited to calculating the state evolution in the Schrodinger 
picture; moreover, it is straightforward to generalize our method to multi-component BEGs 
and to higher-order corrections. 

PACS numbers: 03.75.Be, 67.85.Hj, 03.75.Mn 


I. INTRODUCTION 

We consider the ground state and dynamics of a dilute-gas BEC of N bosonic atoms trapped in 
an arbitrary external potential. In order to describe how interparticle correlations modify the Gross- 
Pitaevskii (GP) equation, we go to the next level of approximation, the Bogoliubov approximation. 
The Bogoliubov approximation [1-3] is important for several reasons: (i) it tells when the Gross- 
Pitaevskii (mean-held) approach begins to break down; (ii) it describes small deviations from the 
Gross-Pitaevskii equation and can be used to study the stability of a BEG; (iii) it enables the 
calculation of how impurities change the behavior of a BEC; and (iv) it is useful for studying phase 
coherence between BECs. 

Conventionally, in the Bogoliubov approximation, the condensate is treated as being close to 
a state in which all the bosons occupy a coherent state of a particular condensate mode (i.e., a 
particular single-particle state). When particle loss is negligible, however, the real condensate is 
much closer to a number state than to a coherent state (see Fig. 1). 

Since a coherent state has a well-dehned phase, the conventional Bogoliubov approximation 
breaks the 17(1) symmetry possessed by the condensate; consequently, a fictitious Goldstone 
mode [4, 5] is present in the Bogoliubov Hamiltonian. Because there is no restoring force on 
the Goldstone mode, the Bogoliubov ground state is not well defined; worse, the Goldstone mode 
causes the condensate state to deviate linearly in time from a single condensate in a coherent state 
(i.e., this is a secular deviation, not an oscillation). This problem is particularly pesky when the 
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condensate is in a trapping potential, where the Goldstone mode is a mixture of the condensate 
mode and modes orthogonal to it and thus cannot be removed easily. The solution to getting rid 
of the unphysical Goldstone mode is to adhere to the fact that the condensate has a fixed number 
of particles, i.e., by using a Bogoliubov approximation where particle number is conserved. 

Many authors have considered the number-conserving Bogoliubov approximation. Girardeau 
and Arnowitt [6] were the first to propose a theory for the ground state and excited states of 
many bosons based on a particle-number-conserving (A^-conserving) formulation of the Bogoliubov 
quasiparticles. G. W. Gardiner [7] introduced a somewhat similar approach to Girardeau and 
Arnowitt’s, but with an emphasis on the time-dependent case; C. W. Gardiner et al. [8] then 
applied this approach to the kinetics of a BEG in a trap. Castin and Dum [9-11] gave a modihed 
form of the Bogoliubov Hamiltonian where the terms that break the ?7(1) symmetry are removed 
by a projection operator. Sprensen [12] generalized the Castin-Dum result to the two-component 
case. S. A. Gardiner et al. [13-15] improved the Gastin-Dum result in the multi-component case 
using an expansion in powers of the ratio of noncondensate to condensate particle numbers. Several 
authors [16, 17] discussed the truncated Wigner approximation, which provides a way to implement 
a number-conserving Bogoliubov approximation in a phase-space description. 

An independent approach is founded on a number-conserving BCS-like ansatz introduced by 
Leggett [18, 19]. Leggett’s ansatz uses the state, 

/if t t 

I^Legg) ot ( “o«o + 2 2^ Afca],a]_fcj |vac) , (1.1) 

k>0 


as a model for analyzing the properties of the ground state of a homogeneous BEG; here A;;. < 1 
and a|, creates a boson of momentum k. Dziarmaga and Sacha [20] generalized Leggett’s ansatz 
to the inhomogeneous case while retaining a similar pair-correlated form. 
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Here Oq is the creation operator for the condensate mode and the creation operators and the real 
numbers are derived from the Bogoliobov Hamiltonian by using a singular-value decomposition. 
Later, Dziarmaga and Sacha generalized their results to the time-dependent case [21], where they 
showed that if the system starts in a Bogoliubov vacuum state, it remains in a state of the same 
structure. The pair-correlated-state approach that Dziarmaga and Sacha introduced is closely 
related to the extended catalytic state approach discussed in the current paper. In [22], we study 
pair-correlated states of the form (1.2) in great detail and demonstrate this equivalence for the case 
where the coefficients A^ are considerably less than one—i.e., the regime where there is a single 
dominant condensate wavefunction; in addition, we derive analytical expressions for the physical 
quantities (particularly the single- and two-particle reduced density matrices) associated with pair- 
correlated states in the large-limit when 1 — A[„ ~ 1/A^, a regime where more than one mode 
can be macroscopically occupied. 

A number-conserving Bogoliubov approximation yields qualitatively different results from one 
that fails to conserve particle number; among these differences are the following. Villain et al. [23] 
removed the “zero-momentum mode” from the Boboliubov Hamiltonian and thereby showed that 
the collapse time of the phase of a BEG is relatively short and, in some cases, vanishes in the limit 
of a large number of atoms. Danshita et al. [24] investigated collective excitations of BECs in a 
box-shaped double-well trap using the number-conserving Bogoliubov approximation. Trimborn et 
al. [25, 26] investigated the artificial number fluctuations in methods that ignore the fixed particle 
number and showed that these lead to ambiguities and large deviations in the Bose-Hubbard 
model. Oles et al. [27] predicted large density fluctuations in a two-component BEG close to the 
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FIG. 1. Phase-space representations for (a) a coherent state with complex amplitude a and (b) a number 
state with particle number N. A number state is distributed phase-symmetrically on the phase space, and 
no definite phase can be attributed to it; in contrast, a coherent state has a well-defined phase. 


phase-separation regime using an ansatz in which the number of atoms in each component is fixed. 
Schachenmayer et al. [28] studied the collapse and revival of interference patterns in the momentum 
distribution of atoms in optical lattices using a number-projection method. Billam et al. [29] 
went beyond the number-conserving Bogoliubov approximation and studied large depletion of the 
condensate by considering the coupled dynamics of the condensate and noncondensate fractions. 

We return to the number-conserving Bogoliubov approximation in this paper and develop a 
particularly transparent method of deriving the relevant equations. Our approach to a number- 
conserving Bogoliubov approximation is to “encode” the many-body wavefunction of the BEC 
in the Wparticle sector of a state we call an extended catalytic state (ECS), by which we mean 
a coherent state for the condensate mode and a state to be determined by the dynamics for 
the orthogonal modes of the atoms. Using a time-dependent interaction picture, we move the 
coherent state to the vacuum, thus making all the field operators formally small compared to 
The resulting Hamiltonian can then be organized by powers of Requiring the terms 

of order to vanish in the interaction-picture evolution equation gives the GP equation for 

the condensate wavefunction. Going to the next order in the evolution equation, N^, we derive 
equations equivalent to those found by Castin and Dum [10] for a number-conserving Bogoliubov 
approximation. In contrast to other approaches, ours allows one to calculate the state evolution 
in the Schrodinger picture, and it also has advantages in considering higher-order corrections and 
extensions to multi-component cases. 

In Sec. II we introduce the ECS (Sec. II A) and the interaction picture in which the conden¬ 
sate mode is displaced to vacuum (Sec. IIB). We then derive the equations that govern the ECS 
dynamics in this interaction picture (Sec. IIC). The GP equation and the Bogoliubov Hamilto¬ 
nian arise naturally as we organize the interaction picture Hamiltonian by powers of We 

make the Bogoliubov Hamiltonian number conserving by adding to it an auxiliary term that does 
not change the dynamics in the Wparticle sector. Finally, we transform back to the Schrodinger 
picture and find a particularly simple form for the ECS dynamics at Bogliubov order (Sec. HD). 
In Sec. HI we generalize our approach to two-component BECs (Sec. HI A) and show how spin 
squeezing is generated in two-component systems (Sec. HIB). Section IV summarizes the results 
of the paper. This paper is based on Z. Jiang’s PhD dissertation (Chapter 3) at the University of 
New Mexico [30]. 
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II. EXTENDED CATALYTIC STATE AND NUMBER-CONSERVING BOGOLIUBOV 

APPROXIMATION 

Quantum optics teaches that coherent states are easier to deal with than number states, and 


what is true there is true here as well. Indeed, the usual mean-field approximation to BEC evolution 
is based on the assumption that the BEC is in a coherent state of a condensate mode [31]. A 
problem with this approach is that the number of particles in a BEC is usually fixed, whereas 


coherent states are superpositions of states with different numbers of particles. A related problem 
is that assigning a coherent state to a BEC breaks its phase symmetry, thus causing problems in 
developing the Bogoliubov approximation. 

A. Encoding the state of a BEC in an extended catalytic state 

Our philosophy for dealing with these problems in a BEC that has a hxed particle number N 
is to extend the BEC state IV’tv) to a state jV’ecs)) which the condensate mode is in a coherent 
state, but the V-particle sector is the same as ji/’jv) within a normalization constant. Consider an 
arbitrary state |'07v) with N particles, for which we have the relative-state decomposition in the 
number basis of the condensate mode. 


N 


I'V’at) — ^ 1^ “ ^)o ® I ^m)± ) ■^±\^m)± — \^m)± 1 (2-1) 


( 2 . 1 ) 


M=0 


where the kets labeled by 0 and _L apply to the condensate mode and to all the modes orthog¬ 
onal to the condensate mode, respectively. The operator A/[l is the particle-number operator for 
the orthogonal modes. The state |C^) for the orthogonal modes, which has M particles in the 
orthogonal modes, is not necessarily normalized. The key to our approach is that the state (2.1) 
can be written as 




( 2 . 2 ) 


where "P^v is the projection operator onto the V-particle sector and 



(2.3) 


is an (unnormalized) state of the modes orthogonal to the condensate mode. 
We now introduce the extended catalytic state, 


l^ecs) = l«)o® |i^)± , 


(2.4) 


which is related to the physical state by 



The extended catalytic state is a direct product of a coherent state |a)o in the condensate mode 
and an unnormalized state |f2)_|_ of the orthogonal modes. Notice that once o. is specified, the 
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extended catalytic state has a one-to-one correspondence with the physical state. The structure 
of the extended catalytic state allows us to study the dynamics of a BEC in the Schrodinger 
picture, and we will see that the structure is preserved throughout the evolution in the Bogliubov 
approximation. 

For a pure condensate with no depletion of the condensate mode, the modes orthogonal to the 
condensate mode are in vacuum, and the overall state has the form 

I'^Tv) = l^)o«'|vac)_L . (2.6) 

In this case we have 

IV’ecs) = l«>o«> |vac)^ . (2.7) 

Generally one expects that a dilute-gas BEC has a state close to that of a pure condensate, in which 
case the noncondensate state |f2)_|_ is close to the vacuum; we want to develop an approximate 
description based on this expectation. To do so, notice that the encoding into an extended catalytic 
state works for any value of a. In other words, one has the freedom to choose a at will; after the 
projection, all values of a yield the same physical state. Nonetheless, we stick to the choice 
\a\ = , for the reason that we make approximations in deriving the dynamics of I’j/’ecs) the 

projection onto the A^-particle sector can amplify the errors due to these approximations. To keep 
these errors under control, we center the number distribution of the coherent state at the actual 
atomic number N. The phase of a is yet another matter, which we discuss further below. 

The BEC Hamiltonian conserves particle number and thus commutes with the particle-number 
operator. As a consequence, the evolution operator hl{t) commutes with 'P^y, allowing us to move 
the evolution operator through the projection onto the Al-particle sector so that it acts directly on 
the extended catalytic state: 

\^N{t)) =^(i)l^iv(0)) 7^jvl^ecs(i)) • (2-8) 

a a 

To find |■J/’^^(^)), one solves for |'0ecs(^)) = ^(^) l'0esc(O)) and then projects onto the A'-particle 
sector. 


B. Interaction picture 

The first step in developing the Boboliubov approximation is to go to an interaction picture in 
which the condensate mode is displaced from a coherent state to vacuum. To do this, we start with 
a condensate mode defined by a time-dependent single-particle state |(/>(t)), which has wavefunction 

(/)(x, t) = (X I (/)(t)) . (2.9) 

The Schrodinger-picture field operator 'k(x) satisfies the commutation relation 

[T(x), 'I'^(x^)] = (I(x — x') . (2.10) 

The annihilation operator for the condensate mode is related to the Schrodinger-picture field op¬ 
erator by 




/ 


0*(x,t)T(x) 


dx = ((/)(t) I T ) 




( 2 . 11 ) 
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Here, in the final two equalities, we introduce a shorthand notation for the integral as bra-ket 
inner products between a single-particle state and the field operator. The creation operator for the 
condensate mode is 

4h) = j ^'^(x)(/'(x,t)dx = (^ I I . (2.12) 

Here and throughout this section, complex conjugation in the single-particle Hilbert space is defined 
relative to the position representation. The bra-ket notation introduced here, though ad hoc, is 
useful for manipulating the complicated expressions that arise as we proceed, more so once we get 
to the two-component case in Sec. HI. Notice that the annihilation and creation operators have 
two different bra-ket forms, both of which are used in our treatment. 

The field operator can be written as 

^(x) = a0(t)</'(x,t)-h ^'_L(x,t) , (2.13) 

where (x, t) is the field operator with the condensate mode excluded. In the Schrodinger picture, 
'k(x) is time independent, but the split between a condensate mode and orthogonal modes 

introduces time dependence because the condensate mode is changing in time; hence, both 
and are explicitly time-dependent operators in the Schrodinger picture. In terms of our 

shorthand notation, we can write 

ivk^(t)) = |T) - = 1 ^) - imm) 1^1^) = Q(t)i'i^), ( 2 . 14 ) 

where 

Q{t) = 1 - P{t) = 1 - \4>{t)){4>{t)\ (2.15) 

is the projector onto the single-particle space orthogonal to the condensate mode, with P(t) = 
\(l){t)){(j){t) \ being the projector onto \4>{t))- Notice that Q*{t) = 1 — P*{t) = 1 — \4>*{t)){(j)*{t)\. 

The extended catalytic state for a pure condensate in the time-dependent condensate mode 
\(t){t)) is 

T» (a, </)(!)) I vac) = |a,())(t))o ® |vac)_i_ , (2-16) 

where the displacement operator P(a, ()>(t)) for the condensate mode, which we usually abbreviate 
as T’(t), is defined as 

V{a, 4>{t)) = V{t) = exp(a4(t) - a a^(t)) , (2-17) 

The state (2.16), which describes a pure condenstate with no depletion, is the one we perturb 
about in developing our approximate description. 

We can now introduce the desired interaction picture as the one where the condensate mode is 
displaced to vacuum; i.e., states transform to 

IV’int(i)) = ^^(a><('(i))ll^ecs(^)) =^int(i) l^int(O)) > ( 2 - 18 ) 

where 

uut) = («, m) m m) ( 2 . 19 ) 

is the evolution operator in the interaction picture. The Schrodinger-picture evolution operator 
U{t) obeys the Schrodinger equation 

m 


dt 




( 2 . 20 ) 
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where T-L{t) is the (possibly time-dependent) BEC Hamiltonian. The time dependence of the 
condensate wavefunction which enters into the displacement operator through 

the annihilation and creation operators, and is to be determined. 

The interaction-picture evolution operator obeys the equation 


ih — ^ = ihV\t)U{t) 11 (0) -|- ihv\t) 71(0) 

CL L {Jib 

= (^ihv\t)V{t) + v\t)'H{t)V{t))u,^,{t) . 

The time derivative of the displacement operator is 


( 2 . 21 ) 


Putting this expression into Eq. (2.21), we have 


( 2 . 22 ) 


a =«!».(*) 

where the interaction-picture Hamiltonian reads 


(2.23) 


'^int(i) = ( \a\^ I ^{t)) + “0(o) +'b>\t)n{t)V{t) . ( 2 . 24 ) 

Equivalently, we have 

= 'Wint(i) I'V’intW) • (2-25) 

In the interaction picture the field operator takes the form 

21^(f) 'I'(x) T>{t) = 'I'(x) -|- a ())(x, t) . (2.26) 

An expansion of 'Hi-ati't) in powers of l/|a| = 1/A'^'^^ is a good approximation as long as the field 
operator T(x) is small relative to the interaction-picture displacement a(f){x,t), i.e., more formally, 
as long as the one-particle density matrix is small in the sense that 

/3int(x,x') = (-0i„t|T'^(x')^'(x)|^int) ~ • (2-27) 


This requirement is satisfied as long as the system is a condensate. We now turn to using the 
expansion in powers of 1/A"^'^^ to derive the number-conserving Bogoliubov approximation. 


C. Number-conserving Bogoliubov approximation 


In second-quantized form, the model Hamiltonian for the BEC is 


m = I 


tVx)( -—V^ + V{x,t) )T(x) -b -4'Vx)4'Tx)^r(x)T(x) 
2m 


h 2 


£ iT/t ,Ttt / 


dx , 


(2.28) 


where the first term is the second-quantized Hamiltonian for particles trapped in a potential V (x, t) 
and the second term represents the two-body scattering energy. The only explicit time dependence 
in the Hamiltonian (2.28) comes from a possible time dependence in the trapping potential V (x, t). 
The present approach is also valid if the interaction strength g is time dependent. Such time 
modulation can be achieved when the s-wave scattering length is controlled by means of, for 
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example, a Feshbach resonance. For our expansion in powers of 1 /|q!| = l/N to work, we must 
have that g\oi'f‘ is of order . 

Going to the interaction picture, we have 


^int W = [ |a| ( 4>{t) I 4>{t )) + a «0(d - « Hit) ) + ^ 


(2.29) 

(2.30) 

(2.31) 
dx, (2.32) 


where we neglect terms of order or smaller. The c-number term (2.30), of order N, is, in the 

time-independent case, the mean-field energy of the BEC; its only effect, in general, is to introduce 
a global phase, so we ignore it henceforth. 

By requiring the linear term (2.31), of order to vanish, we get 

/ \ 

ih(j){x,t) = (- — V‘^ + V{x,t) + g\af\(j){x,t)f j (j){x,t) = Hgp{t)(j){x,t) , (2.33) 

which is the celebrated Gross-Pitaevskii equation. The single-particle GP Hamiltonian is 

+ V{t) + g\af . (2.34) 

The structure of our approach is now clear. By going to the interaction picture, the mean-field, 
Gross-Pitaevskii evolution is removed, and then by neglecting the terms of higher order than 
we are left with the quadratic Bogoliubov Hamiltonian 


qj — 
/Lbog 




gp 




dx 


(2.35) 


To display the symplectic structure of the Bogoliubov Hamiltonian, we write it in the matrix 
form 



((T| (4/t|)Hbog 



(2.36) 


where the colons denote normal ordering of annihilation and creation operators, and the 2x2 
matrix reads 


^ /H^^ + g\a\^\cP\^ ga^ \ 

^ V ^^gp + 5|a|VlV 

Notice that the normal ordering has an effect only on the lower-right corner of the matrix Hbog- 
As shown by Lewenstein and You [32], has a nilpotent subspace, where phase diffusion 

takes place. Such phase diffusion is not physical, but rather is a consequence of the arbitrary phase 
assigned to the condensate wavefunction, i.e., to a. This problem was addressed by introducing 
number-conserving approaches [6, 7, 10]. Particularly in the work of Castin and Dum, a systematic 
expansion of the field operators was used in deriving the equations for the number-conserving 
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Bogoliubov approximation. The aim is to eliminate the artificial nilpotent subspace that gives rise 
to the phase diffusion. Here we solve the same problem by introducing an additional contribution 
to the Hamiltonian, an auxiliary, explicitly time-dependent, Schrodinger-picture Hamiltonian 
which does not affect the A^-particle sector of |■0ecs(^)) thus keeps the physical state |'07v(^)) 
unchanged, i.e., 

iP^J-(t)|^ecsW> =0. (2.38) 

With this term T'(t), we can solve the phase diffusion problem by eliminating the nilpotent subspace 
of b/j^Qg. 

To determine the form of T'(t), we must go to the Bogoliubov level of approximation, but for 
now let us suppose T'(t) takes the form 

Ht) = -'^{Af-Nf + (a4(i) - N) - N) . (2.39) 

Here 

AA = / ^^(x)^(x)(ix = + j t|j_(x, t)^'_L(x, t) dx (2.40) 

is the total particle-number operator, and A/j_ = M — is the particle-number operator for 

all the modes orthogonal to the condensate mode, i.e., the depletion number operator. The time- 
dependent parameter r]{t), which is to be determined, is of order N~^. The operator j^j_, also 
to be determined, is of the order and is a linear function of the annihilation and creation 

operators of the modes orthogonal to the condensate mode; it thus commutes with and a^. 
The first term in Eq. (2.39) clearly satisfies Eq. (2.38). For the other two terms, we have 

0 = T’tv (AA — ~ 'PN T 1. ~ |q^)o®I^)±) (^•‘^1) 

and 

0 = (|a|2_iV) = {a* - N)Fi\a)o®\Q) i. , (2.42) 

where in the first equation we use aaj, -|-A/j_ — N = 00 ( 0 ; — a^) +N — N. As long as the condensate 
mode stays in a coherent state with amplitude a, these two terms do not affect the physical state 
|'0jv(t)). We show that the condensate mode does remain in a coherent state at Bogoliubov order 
in Sec. H D. 

An astute reader will have noticed that the auxiliary Hamiltonian (2.39) is not Hermitian. This 
is not a problem at Bogoliubov order, however, because the only nonHermitian term in F{t) is 
A/j_Tj_, which, being of order , can be neglected in the Bogoliubov approximation (order N^). 

Going now to the interaction picture, we have 

-^int(i) =P\oi,4>{t)) F{t)V{a,cl){t)) = -^{aa\ + a + {aa\+M^)F^ +a a^F\ , 

(2.43) 

where the identity \a\ = N is used to cancel several terms. If we now discard terms of order N ' 
or smaller (in doing so, recall that interaction-picture field operators are order N^), we obtain 

-T— ^/r)||2j' I 2't'j'|/*\2 f .7— I =(= ^ I 12 / o /I /I ^ 

*^int ^ ^ “1“ ) *^0^0 j “1“ *^J_ ^ ^0 *^J_ ^ ^ (2.44j 

+ a{-^\c^)F^ + a{(t)\-^)F{-'^\a\^ . 


(2.45) 
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Here we normally order the creation and annihilation operators of the condensate mode in prepara¬ 
tion for incorporating into the main Bogoliubov Hamiltonian; this normal ordering introduces 
the c-number term —•r]\a^j‘l. This term could be important as a second-order correction to the 
condensate energy, but it only adds an overall phase to the evolving quantum state, so we neglect 
it henceforth. In Eq. (2.45), we introduce the bra-ket notation of Eqs. (2.11) and (2.12). The 
modified (number-conserving) Bogoliubov Hamiltonian then takes the form 


"^ncb = ^bog + -^int ■ (2.46) 

To eliminate the phase diffusion, we choose 

^(i) = 9 j= g{4>\\(p\^\4>) = g{(j)*\\(j)\^\(j)*) = g{(p*= g{4>\4>'^\(f)*) (2.47) 

and a Hermitian 

Tj_{t) = -ga* J gaj ^(x, f)|(/)(x, f)|^T|j_(x, f) dx (2.48) 

= -ga*{(j)\\(f)\‘^Q\^) - ga{(j)\(j)‘^Q*\^'^) = -ga*\Q*- ga{^\Q\(i)\‘^\(j)) . (2.49) 


It is now a tedious calculation to show that 

T- _ ^ ol |2 t I 2 t t , \ 

‘^int 2 V ^ J 

- g(^a\^a\ +{af j (j)*{x,t)\(j){x,t)\^'ii{x) dx + R.c^ 

= -| (^2|a|^ {^\P\(j)\'^P\^) + {a*f {'^^P*{(f)*fP\^) +a‘^ 
-g(\a\‘^ (^|P|</.|2q|4') + 

+ (af {^^\Q*{(j)*fp\^) + a {^\P(j)^Q*\^^)^ . 

Translating this into matrix notation, we get 

-^int = 2 : ((^1 (^^l) ^int • ’ 

where the 2x2 matrix is 

/ \a\\Q\cl^\^Q-\<P\^) \ 


(2.50) 


(2.51) 


(2.52) 


(2.53) 


In Eqs. (2.47), (2.49), and (2.51), we use the bra-ket notation, which is the easiest way to carry 
out the algebraic manipulations; for this purpose, it is useful to notice that = 

('kI I ) and its conjugate, = ('k^|P*(())*)^|'k). The bra-ket manipulations 

generalize straightforwardly to the two-component case considered in Sec. HI. 

The number-conserving Bogoliubov Hamiltonian (2.46) now reads 


n 


neb 


1 

2 



Q ■ 


• ^ 


(2.54) 
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with 


L-l =1-1 -1- F = 

^ncb ^bog ' ^int 


+ g\a\‘^Q\(f>\‘^Q 


ga^ Q4>^Q* 


gp 

g{a*fQ*icl>*fQ + g\a\^Q*\cl>\^Q* 


(2.55) 


This number-conserving Hamiltonian is the same as that found by Castin and Dum [10] using a 
systematic expansion of the field operators. The two approaches give the same dynamics for the 
Bogoliubov approximation and thus are equivalent to order The difference between the two 
approaches is that Castin and Dum derive their results in the Heisenberg picture, whereas we use 
the Schrodinger picture and a closely related interaction picture. To summarize our approach, we 
move the coherent state of the condensate mode to vacuum by going to a time-dependent interaction 
picture. The interaction-picture Hamiltonian, when organized by powers of iV , gives the GP 
equation at order and the conventional Bogoliubov Hamiltonian P^bog order . The 

conventional Bogoliubov Hamiltonian suffers from the artificial problem of phase diffusion because 
a zero-momentum mode arises from choosing a phase for the condensate wavefunction. Without 
affecting the Wparticle sector, we remove the phase diffusion problem and derive the number- 
conserving Bogoliubov Hamiltonian (2.55) by introducing the auxiliary term (2.39). 

It is useful below to divide the number-conserving Bogoliubov Hamiltonian (2.54) into its two 
natural parts, = ^gp + Here 


^gp 2 • (^^1 (^i^t J • 

= (^l^gpl^) = «W('^l-^gpl'?^) + 4('^l-^gpl^-L) + (^±l^gpl<?^)«0 + (^ul^gpl^u) > 

with 


H 


gp 


///gp 0 \ 

V 0 hJ ’ 


is the GP part of the number-conserving Bogoliubov Hamiltonian, and 


/C = 


1 

_ • 

2 • 






with 

^ / \a\^Q\4>\^Q aQ(t?Q* \ ^ 


(2.57) 


(2.58) 


(2.59) 


describes the additional coupling of the orthogonal modes coming from two-body scattering. Notice 
that once 1C is written in the second form of Eq. (2.58), we can omit the projectors Q and Q* from K. 

We find it useful to introduce an orthonormal basis of single-particle states at t = 0, {|Xj(0))}. 
We choose |xo(0)) = |'/(0)); the j > 1 states are then a complete set of modes orthogonal to |(/(0)). 
We evolve all these states forward in time using the GP Hamiltonian Hgp{t), i.e., 

\Xj{t)) = Ugp{t)\xji0)) , (2.60) 


where the single-particle GP evolution operator Ugp{t) is the solution of 


ih 


dU^p{t) 


dt 


-Hgp(t)Hgp(t) . 


(2.61) 
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The condensate mode |<^(t)) = |xo(^)) satisfies Eq. (2.60) by virtue of the GP equation (2.33). 


Notice that Ugp{t)Q{0)Ugp{t) = Q{t). 

The corresponding annihilation and creation operators are 

SW = (v(<)l'l'> = ('l'’iq(()>, (2.62) 

»;(i) = {xK*)l'I'')=('I'lv(t)>. (2.63) 

The field operators with the condensate mode excluded can be written in our bra-ket notation as 

\'^±(.t)) = '^aj{t)\xj{t)) , (T|L(t)| = ^aj(t)(x*(t)| , (2.64) 

i>i i>i 

j>i i>i 


In terms of these time-dependent single-particle states, the coupling Hamiltonian (2.58), with all 
the time dependence indicated explicitly, takes the form 

K:{t)=g\af «i(^)«fcW(Xi(i)|l0l^(i)|Xfc(i)) 
j,k>l 

j,k>l 

( 2 . 66 ) 


D. Dynamics in the Bogoliubov approximation 

We turn now to the dynamics of the BEC within the Bogoliubov approximation. We begin 
by recalling that in the Bogoliubov approximation, an approximate interaction-picture evolution 
operator is constructed from the number-conserving Bogoliubov Hamiltonian (2.54), i.e., 

ifi = ^ncb(i) ^int(i) , (2-67) 

The corresponding interaction-picture evolution is |'0int(^)) =^int(^) l'0int(O))' 

Our first task is to confirm that if the condensate mode begins in a coherent state, it remains in 
a coherent state with the same complex amplitude a under the evolution of the number-conserving 
Bogoliubov Hamiltonian (2.54). Eor this purpose, it is instructive to use Eqs. (2.56) and (2.58) to 
divide the field operators in the number-conserving Bogoliubov Hamiltonian into a contribution 
from the condensate mode and a contribution from the orthogonal modes: 

'^ncbii) = + al(t^\ffgp\'^±} + (^'±|i?gp|0)o</, + Kcbx _ 

( 2 . 68 ) 

= a^(())|77gp|'k) -I- {'if±\Hgp\4>)a^ -I- 'Hnchx ■ 

Here 

Kcb±(t) = I : (('I^ul ('kii)Hbog : = ('k^l^gpl^u) +/C (2.69) 

is the Hamiltonian for the orthogonal modes; we can use instead of Hn^b because the projectors 
Q and Q* have no effect. Using the GP equation (2.33), we can rewrite Eq. (2.68) as 

"^ncb T f^{ ^_L I 0 T "^ncb-L ? 


(2.70) 



13 


from which we can immediately verify the commutator identity 

In the Heisenberg picture, we have 

([«</-(*)> ^ncb(^)] + = 0 , (2.72) 

from which we conclude that 


^iLt(0«<^(t)^int(i) = a^(0) . (2.73) 

The conservation of in the Heisenberg picture implies that in the interaction picture, if the 
condensate mode begins in vacuum, it remains in vacuum. Equivalent to this statement is the 
statement that in the Schrodinger picture, the condensate mode is always in the coherent state 
^(a, i?i)(t)) |vac)o = |a,(/>(t))o- As a result, Eq. (2.38) is always satisfied, and the auxiliary Hamil¬ 
tonian T" does not affect the physical state. Notice that this means that in the extended catalytic 
state, the condensate mode does not become entangled with the other modes. When we project the 
extended catalytic state to the Wparticle sector to obtain the physical state of the BEC, however, 
entanglement makes its appearance. 

To get into the Schrodinger picture, as we promised to do, requires some additional formal 
apparatus. The effort is worthwhile, however, because it reveals the role of the GP part of the 
evolution at Bogoliubov order and identifies the most instructive formulation of the Bogoliobov 
evolution. We begin by defining an evolution operator 7/gp(t), which changes according to the GP 
part of the number-conserving Bogoliubov Hamiltonian, 

= (^|^gp(t)|'k)^7gp(t) . (2.74) 

It can be seen that 

^^^^ip(i)^(x)7/gp(t) =7/^p(t)[T(x), 77gp(t)]7/gp(t) = j dx'i7gp(x,x',t)7/^p(t)T(x')7/gp(t) , 

(2.75) 

whose solution can be written in terms of the single-particle GP evolution operator Ugp{t) [see 
Eq. (2.61)], 

7/^p(t)4'(x)7/gp(t) = j dx't/gp(x,x',t)T(x') ; (2.76) 

in more symbolic form, we have 

Z7gtp(t)|T)^gp(t) = t/gp(t)|T). (2.77) 

This is the unsurprising conclusion that the field operator evolves under the GP part of Bogoliubov 
Hamiltonian according to the single-particle GP evolution. Notice that 7/gp(t) does not mix creation 
and annihilation operators. 

Consequences of Eq. (2.77) are the following: 

z^gV^)«0w^gpW=^ipW(<^w I 'J^)^gpW = imugpitm = {m i 'k) = a^(o), ( 2 . 78 ) 

z7gVt)|T^(t))z7gp(t) = Q{t)uiitmu^p{t) = = u^pim±m . 

(2.79) 
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Equation (2.78) says that is conserved under the GP evolution. More generally, Eq. (2.79) is 
the statement that the annihilation operators for all the orthogonal modes propagated using the 
single-particle GP Hamiltonian [see (2.62)] are also conserved, just as in Eq. (2.78), i.e., 

Ulp{t)aj{t)l(gp{t) = aj{0) = ttj . (2.80) 

The evolution operator U^p{t) is very close to being a formal device: it translates between the 
natural modal descriptions that apply at different times as the single-particle states evolve under 
the single-particle GP Hamiltonian. Indeed, it is easy to see that 

^gpWho></>(0);ni,Xi(0);n2,X2(0); ■ • • ) = 1%,'/>(*); Xi(i); X2(i); • • • ) , (2-81) 

where \nj,Xj{t)) = [aj(f)]"''’'|vac)/yTi)^ is the state with rij particles in the single-particle state 
\Xjit))- This gives us U^{t) as an explicit basis transformation: 

^gpW = X] ho></>W;«iiXi(i);?^2,X2(i);■ ■ • )(’^o;</'(o);«iiXi(o);n2,x2(o);..■ I ■ (2.82) 

no,ni,n2,... 

We hnd it useful to have available the restriction of 7/gp(t) to the orthogonal modes: 

^gp±W = X] hiWi(i);?^ 2 ,X 2 W; • ■ ■ )(f^i,Xi(0);n2,X2(0); ■ ■ • I ■ (2.83) 

ni.nj,... 

All this suggests going to an interaction picture relative to the GP part of the Hamiltonian and 
solving for the evolution operator 


which obeys the evolution equation 


where 




= nmt) 


ic{t)=ul{t)icm^^{t) = i: ((^±(0)1 ('h 


m 


|T^(0))' 

l'I^i(O)), 


Here the matrix of symplectic structure is 

0 

^gp 

'q(o) 0 

, 0 Q*(0) 






VO ^gpW. 


= 9 


'uUt) 


0 


0 uUt) 




0 


0 


'Q(O) 0 ^ 

, 0 Q*(0)^ 


(2.84) 

(2.85) 

( 2 . 86 ) 


(2.87) 


The projectors <3(0) and Q*(0) can be omitted when the matrix K{t) is inserted into Eq. (2.86). 
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We can get a better idea of what the Hamiltonian (2.86) means by writing it in terms of the 
time-dependent single-particle states considered in Eqs. (2.60), (2.62), and (2.63): 

iC{t)=g\a\^ ^ a]ak{xj{0)\ulp{t)\(J)\‘^{t)Ugp{t)\xk{0)) 
j,k>i 

+ f Z] ((«*)^«i«fc(Xi(0)|C/J,(t)(</>*)^(t)t/gp(t)|Xfc(0)) (2.88) 

j,k>i 

+ a‘^a]al{xj{0)\ulp{t)(p‘^{t)U*p{t)\xli0))'^ . 

= 9\oif a]ak{xjit)\\(p\‘^{t)\xkit)) 

(2.89) 

j,k>l 


The form (2.89) can be obtained directly from applying Eq. (2.80) to Eq. (2.66), or it can be 
obtained by the route through Eq. (2.88), which shows that the role of the single-particle GP 
evolution operators is to transform the coupling matrix elements into the time-dependent basis 

{\xjm- 

We can now write the Schrodinger-picture evolution operator as 

^/(t) =P(a,<()(t))^^i„t(i)^?^(«,</-(0)) =2?(«,</>(t))^gp(t)V(t)pt(«,</.(0)) . (2.90) 

Using Eq. (2.78), we have h(gp{t)'D(^a, 4>{t))h(gp(t) = Vi^a, , and noting that V(t) only acts 
on the orthogonal modes, we can remove the displacement operators from the evolution operator, 
obtaining 


l({t)=U^p{t)V{t). (2.91) 

The upshot of all this is that the Schrodinger-picture evolution involves, first, evolution of the 
orthogonal modes, with fixed creation and annihilation operators, under the Hamiltonian (2.89) 
and, second, translation of the mode structure to the time-dependent modes evolved using the 
single-particle GP Hamiltonian. The displacement of the condensate mode to vacuum is, as we 
anticipated, a formal device for developing the expansion in powers of l/|a| = it disappears 

from the final Schrodinger-picture evolution. 

Suppose now that, in accordance with our general assumptions, the initial extended catalytic 
state is |^ecs(0)) = |a, (/)(0))o(8) |12(0))_l. Then V(t)| V’ecs(O)) = l«> '/>(0))o® V(t)| 12(0 ))_l, since V{t) 
only acts on the orthogonal modes. The operator hl^p{t) translates this state to the modes that 
apply at time t, giving 


ll^ecs(O) = |a,</>W)o® l^(i))_ 


(2.92) 


where 


\m)±=u^^At)v{t)\m). 


(2.93) 


To find the physical state at time t, one projects the extended catalytic state onto the Wparticle 
sector, as specified by Eq. (2.8), which gives 


N 


IV’iv(i)) = \n - 


M=0 


(2.94) 
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where 


1 / IV! 

(jV - M)\ 


(2.95) 


Here projects onto the M-particle sector of the modes orthogonal to the condensate mode 

at time t. Notice that M' does not have to equal M because V{t) is not number conserving. 


III. TWO-COMPONENT BECS 

In Sec. II we discussed how to derive the number-conserving Bogoliubov approximation for a 
single-component BEC by going to an interaction picture where the condensate mode is displaced 
to vacuum. In this section we show that it is a simple task to generalize our method to multi- 
component BECs. We do the two-component case as an example, but the generalization to many 
components is straightforward. 

We are certainly not the first to consider a number-conserving Bogoliubov approximation for 
the multi-component case. Sprensen [12] generalized the Castin-Dum result to the two-component 
case, and this facilitated discussions on spin squeezing in BECs [33-36]. S. A. Gardiner et al. [13- 
15] improved the Castin-Dum and Sprensen results by using an expansion in powers of the ratio of 
noncondensate to condensate particle numbers, which is advantageous for large depletion. Com¬ 
pared to these previous studies, our approach for the multi-component case is distinguished mainly 
by the ability to carry over the single-component case with very little modification, essentially a 
generalization to a spinor notation for the several components. 


A. Number-conserving Bogoliubov approximation for two-component BECs 


In the two-component case the condensate wavefunction, which is generally a single-particle 
state that is entangled between the translational and internal degrees of freedom, takes the form 

\Ht)) = “ («lWI</>l(O)®|l)+«2(i)l02(i))® |2)) , (3.1) 

a ^^ a V / 


where a, which takes on values 1 and 2 in the two-component case, labels the hyperfine levels and 
where \ai\^ + |a 2 |^ = |a|^ = N, with N being the total number of particles. The states |1) and |2) 
are internal states of the bosonic atoms, which we refer to as hyperfine levels because that would 
be a typical situation in a dilute-gase BEC. In the subspace spanned by (8) |1) = |l,<?!>i) and 
[(^ 2 ) ® 12) = [2, (/) 2 ), the single-particle state that is orthogonal to the condensate mode is 


Notice that 


= At ( Oilit)\Mt)) ® |1) - 

a \ 


l>i(i)> = \Mi)) ® |1) 
2,<('2(i)> = IMt)) ® |2) 



alit) 



(3.2) 

a J 

(3.3) 

a / 

(3.4) 


a 
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The field operator that destroys a particle in internal level a at position x is 'I'^(x). In our 
shorthand bra-ket notation for field operators, we have 

= (x| = (cr,x| T) . (3.5) 

In the final form, we extend our notation by introducing a total field operator 

= ( 3 . 6 ) 

CT 

which is a spinor field operator, including both spatial and internal degrees of freedom. It gives the 
hyperfine-level field operators according to (it| T) = iTg-); notice that since = ('^'o-1 x) = 

( T I fj, X), we also have ( T | a) = |. The spinor representation is 

T(x) = (x| T) = ^^'^(x)|(t) . (3.7) 

CT 

We can also write 'ko-(x) = ( | x) = ( | (T,x) and 'I'ct(x) = I ) = (c^^x | ). 

The annihilation and creation operators that destroy or create a particle in internal level a with 
spatial wavefunction ip{x) are 

ba,i} = j ^a(x) dx = ( V' I ) = ( O', V’ I ^ ) , (3.8) 

= j V'(x) dx = ( I V’) = ( ^ I 0 -, V') • (3.9) 

The annihilation operators for the entangled states \(j)) and \(j)) are thus 

= (</' I ^ I 'I'l) + Q^2( '/’2 I '^'2 )) ~ + ^2^2,(^> 2 ) ’ (3.10) 

«</. = (01 ^ (a2(0i 1 ^ 1 )- ai(02 1 ^ 2 ))=^ («2&i,0i - oiib2,4,^) . (3.11) 

The field operator for the atoms in hyperfine level a can be written as 

^<7(x) = + ^a±(x,t) . (3.12) 

The total field operator can be written in a variety of forms, 

1^) = X] ® 1^) + (3-13) 

(J 

= a^it) 1 4>{t )) + d^(i) 10(t)) + I ^'jL (t)) (3.14) 

= a^it)\4>{t)) + \^±{t)) , (3.15) 

where 

P' 16 ) 

(7 

is the total field operator with modes \(p) and |0) removed and 

l^±(^)) = + l'J'jL(i)) = 1^) - a^{t)\<Pit)) = <3(i)l^) 


(3.17) 
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is the total field operator with only the condensate mode removed. The projectors onto and 
orthogonal to the condensate mode, 

P{t) = m)){cl){t)\, Q{t) = l-P{t), (3.18) 

are defined as in the single-component case [cf. Eq. (2.15)]. By using our bra-ket shorthand, all 
the manipulations for two components can be made identical to that for a single component. 

Just as in the single-component case, we perturb about the extended catalytic state for a pure 
condensate that is in a coherent state for the condensate mode: 


T>(a,(/)(t))|vac) = |a, (/)(t))o 0 |vac)_L , 

The physical state is obtained by projecting onto the A^-particle sector. 
In the two-component case the model Hamiltonian for the N atoms is 

nit) = J2 f + VAt)]'^.d^ + y^hL,,J 


(3.19) 


2 m^ 


cr.r 


2 

a.T 


dx . 


(3.20) 


The diagonal terms of the Hermitian matrix give the energies of the internal levels, and the 
off-diagonal terms give the single-particle coupling between the two levels. The real, symmetric 
matrix describes the scattering of the atoms in each component off one another and the cross¬ 
scattering between components. Since the single-particle terms are trivial to treat, the really new 
effect comes from the cross scattering described by gi 2 - 

The next step is to go to the interaction picture where the condensate mode is displaced to 
vacuum, just as in Eq. (2.18). In this interaction picture, the field operators transform according to 

(a, (/)(t)) T^(x) V{a, (j){t)) = a^{t) 0^(x, t) + T^(x) , (3.21) 

thus allowing an expansion in powers of = l/|a|. We can write this transformation more 

abstractly as 

\v\a,<P{t)) ^V{aA{t))) = am)) + |T) , (3.22) 

The interaction-picture Hamiltonian, as in Eq. (2.24), is given by 

'^int(^) = ( l«l^( I ) + a ('h I m)) - a* {m) I Tdit) • (3.23) 


The time derivative of the condensate state is 


^ ® |i) + IMt))) ® |2)) • 


(3.24) 


Putting all this together, we get the interaction-picture Hamiltonian to Bogoliubov order, i.e., 
order 


'^int(i) = [yarn (^-ih^+H„ + ^ygmTf\(t>T\jam + J2f^aTar(l)^ 


+ / ( + Ha + y9aT\ar\‘^ \(l)r\^)am pyrnrarfl^. 


dt 


dx (3.25) 

-|- H.c.^ dx 
(3.26) 




+ J2gmTf\(I^T\j^a + J2y{f>'^aT+9aTaa(t)aa*r(l)*r)'^-, 

/T X T / rr n- 


+ i^J2(yy9aTamam + H.c. 


dx , 


(3.27) 
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where the single-body translational Hamiltonians are 


H = —+ V 
"" 2m„ 


(3.28) 


Just as for a single component, we can neglect the c-number, mean-field-energy term. By 

1/2 

requiring the term of order N ' = |q;| to vanish, we get a pair of coupled GP equations, 

= 0 . (3.29) 

\ j- / j. 

Notice that these are best thought of as coupled equations for the unnormalized wavefunctions, 
ai(pi and Q!202- If is often convenient to have the two GP equations written out separately as 


'ih— -|- Hgp ^ OiiC^i hiOi2 0!24>2 — 0 , 

d 


(3.30) 


+ -figp ) «2<?i2 + ^21 = 0 , 

where the GP Hamiltonians are 

= Hi-\-hwii + gii\ai\ \(j)i\ +5i2|«2l I'?i2l 
= H2 + hW22 + 522l®2l^l<?i2l^ + 92l\oilf \4>lf 


(3.31) 

(3.32) 


(remember that a; 2 i = W 12 and 521 = 512 )- If is also convenient to make the equations compact by 
writing them in terms of spinors relative to the two hyperfine levels so that we can take advantage 
of our bra-ket notation. 


d 

— iH~p— “h 


dt 


\a 2 cP 2 ) 


= 0 , 


(3.33) 


where 

/ fiu \ 

= L "" (2) = hS\1){1\ + ^£/2)(2| + h(Ui2|l)(2| +ha;2i|2)(l| . (3.34) 

\fuj 21 Hy J 

Recognizing that the spinor in Eq. (3.33) is the spinor representation of the state a\(j)), we can 
write the coupled GP equations in the very compact form 

d 

+ -flgp)l<?i) = 0 , (3.35) 

where it is assumed, as our formalism requires, that a does not change in time. 

The coupled GP equations (3.31) and (3.32) imply that 

^(|ail^(<(>il</>i)) =-^(|a2l^(</>2 I '/>2)) = 2Im(a;i2ala2(</>i ^ 2 )) , (3-36) 

*Il^(ai«2(<(>l I </>2» = «l«2(</>l|(l^gp - ^Igp )|<^2) +Il^]L2(l«ll^(</'l \4>l) - \a2?{(t>2 \ ^ 2 )) ■ 

(3.37) 


2 

The first of these ensures that |a| (</>!</>) is conserved; since we require a to be a constant, we 
have that {(j)\^) is conserved, as is implied directly by the compact GP form (3.35). Moreover, 
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our formalism assumes that | (pi) and | (f> 2 ) remain normalized to unity, implying that any temporal 
changes in \ai\^{(pi \ (pi) and |a 2 p ('^2 I '/’2 ) ^re incorporated into the magnitudes |q!i|^ and |a 2 p; 
this simplifies Eqs. (3.36) and (3.37) to 

= 2Im(a;i2ala2(</>i ^ 2 )) , (3-38) 

ih^{aloi2{(pi \(p2)) = oi*iOi2{(pi\{Hf^ - i7^p^)|<^2) +^‘^i2(|ail^ - \(^2?) ■ (3-39) 

Notice also that we can always move any phase changes in and a 2 into \ pi) and |i?i 2 ); means 
that we can always choose and 02 to be real. 

If the internal levels are eigenstates of the single-particle Hamiltonian, there is no single-particle 
coupling of the internal levels, i.e., W 12 = 0 J 21 = 0. One often uses transient, strong coupling of the 
internal levels to induce transitions between the internal levels. This occurs on timescales much 
shorter than that of the nonlinear terms in the GP equation and can be treated separately as 
a sudden single-particle effect while ignoring the nonlinear terms; the result is a sudden change 
in ai and 02 while i;^i(x) and <p 2 {'^) remain unchanged. Thus the single-particle coupling terms 
can generally be omitted when analyzing BEG dynamics; we retain them for completeness in our 
general development of the Bogoliubov Hamiltonian. 

Before moving on, however, we note that if UJ 12 = W 21 = 0, Eq. (3.38) implies that |q!i| and \a 2 \ 
are constant in time. Since we can move any phase changes in ai and 02 into |(/)i) and |(/> 2 ), we 
can assume that and a 2 are constants, which simplifies Eq. (3.39) to an equation for the change 
in the overlap of \pi) and |<^ 2 )^ 

i^j^{(Pi\(t>2) = -H^gp)\4>2) ■ (3.40) 

Furthermore, it is easy to see from Eq. (3.33) that under these circumstances, \(p{t)) satisfies the 
compact GP equation: 


d 


\p) = 0, if UJ 12 = u}2i = 0. 


(3.41) 


The Bogoliubov Hamiltonian governing the dynamics in the interaction picture is given by 
Eq. (3.27). In 4 X 4 matrix form, we have 


^hog=\: (4'2| (^ll (^2l)^bog 


1^2) 

I'I'l) 

\\H)J 


where the matrix f/^og takes the form 


(3.42) 



^12 + 912(^1^2 4>i4>2 

5ii®i 4^1 

512Q^1®2 4^i4>2 

f^21 + 5210^10^2 Ph'p2 

+ 922\0i2?‘\<p2?‘ 

521«1«2 4>l‘p2 

922<^2 P’l 


* * i* 

5210^1^12 P1P2 

+5ii|aiH<?^il^ fi(^2i + 92ia*ia2p*ip2 

* * 

\ 512«1«2PiP2 

522 {alfifplf 

^12 + 512^102 </>l<?^2 H^p + 922\a2f\4’2\^ 


(3.43) 
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To get back to the compact spinor notation, we introduce, along with the matrix (3.34), two other 
matrices that operate in the spinor space defined by the hyperfine levels |1) and |2): 


1 /0 \ 
"VO a2(/)2/ 


With these matrices, we have 


G 


(9ii 9 i2 
\fl'21 922. 




(3.44) 


(3.45) 


Notice that since $ is diagonal and G is real and symmetric, and are both Hermi- 

tian, and they are transposes and complex conjugates of one another; 4>G<1> and are both 

symmetric, and they are complex conjugates and Hermitian conjugates of one another. Using our 
total field operator and interpreting the 2x2 submatrices as operators in the space of the internal 
levels, we can write the Bogoliubov Hamiltonian in the suggestive form, identical to that for a 
single component, 



((^1 (^'^l)^bog 



(3.46) 


To eliminate phase diffusion in the condensate mode, we now introduce the auxiliary (nonHer- 
mitian) Hamiltonian T" in exactly the same form it has in the single-component case [cf. Eq. (2.39)], 

Ht) = -^ (W - iV)2 + (««$,(*) + A4(t) - - iV)4 , (3.47) 

where A/j_ = M — The coefficient rj and the operator Jyi = defined in analogy to the 

single-component case, 

/ 9aT\aafWr\^\(t>a\^\(l>T?dyi (3.48) 

.,r l«l J 

= (0|$G4>*|(^) = = {^\<^G^\cl)*) , (3.49) 


and 


J^± 


\a\ 




—a 


|$G$*Q|T) - = -a 


dx. 


(3.50) 

-a(T|g$G$*|(/>) , 

(3.51) 


where QIT) = |T_|_) is the total field operator with the condensate mode excluded [see Eq. (3.17)]. 
As in the single-component case, t] is of order 1/N and J-j_ is of order The argument that 

the auxiliary Hamiltonian does not change the evolution in the Wparticle sector, as long as 
the condensate mode stays in a coherent state with amplitude a, is the same as that given in the 
single-component case in Sec. HC. 

The transition to the interaction picture goes exactly as in the single-component case, yielding 
Eqs. (2.44) and (2.45) at Bogoliubov order N^. Dropping the c-number term from that result, we 
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find the analog of Eq. (2.51): 


•^int 


- (^|ap (^|P$G^>*Q|^') + |a|^ (^'|Q^>G^>*P|^) 


Translating this result to matrix form of symplectic structure, we have 


-^int(i) = ^ : ((^1 Fintit) 



F-mt — 


where 

\af{Q^G^*Q - $G«>*) {Q<^G^Q* - $G^>) 

^(a*)^(Q*$*G^>*Q - ^*G^*) \af{Q*^*G<^Q* - 
The number-conserving Bogoliubov Hamiltonian matrix assumes the form 

|M/) 


qj _ _ 
^ncb 


^ : ((T| h/ncb 




with 


u — A-/ -\- F = 

^ncb ^bog ' int 


Hgp + \afQ<^G<!>*Q aQ<^G^Q* 

, {afQ*<^>*G^*Q H*p + \a\^Q*^*G^Q* 


(3.52) 


(3.53) 


(3.54) 


(3.55) 


(3.56) 


Using the same strategy as in the single-component case, we have derived the number-conserving 
Bogoliubov Hamiltonian (3.56) that governs the dynamics of a two-component BEC in the inter¬ 
action picture. This Hamiltonian has the same form as the Hamiltonian (2.55) that applies in the 
single-component case; the difference is that here 4>, G, Q, and H^p are themselves matrices. Our 
result conforms with Eq. (3.17) in [12], but in a more compact form. This compactness is a major 
advantage in generalizing to the multi-component case. 

As in the single-component case, it is useful to divide the number-conserving Bogoliubov Hamil¬ 
tonian (2.54) into its two natural parts, Pncb = ^gp + The GP part is 


n 


gp 


1 

2 • 


((T| (T'f|)^gp 



(3.57) 


= (^l^gpl^) = ala^{(f)\Hgp\4>) + a^(</>|Pgp|T_L) -F ±\Hgp\(j))a^ + (T_L|Hgp|4'_L) , 

where 


^gP = 


H. 


gp 


0 h: 


gpy 


The additional coupling of the orthogonal modes, coming from two-body scattering, is 


/C = ^:((T| (Tt|)/< 





(3.58) 


(3.59) 
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where 

\ [q 0 

a\^^*G^) \0 Q* 

(3.60) 

The demonstration that at Bogoliubov order, if the condensate mode begins in a coherent 
state, it remains in a coherent state can be repeated word for word from the single-component case 
considered at the beginning of Sec. IID, and this shows that the auxiliary Hamiltonian J-{t) of 
Eq. (3.47) does not change the evolution in the A^-particle sector. 

The Schrodinger-picture dynamics can be developed using exactly the same approach and sym¬ 
bology we used in the single-component case, beginning with Eq. (2.74) and running through the 
end of Sec. IID. In particular, the Schrodinger-picture evolution operator can be written as in 
Eq. (2.91). The operator Z7gp(t) evolves according to the GP part of the number-conserving Bo¬ 
goliubov Hamiltonian, as in Eq. (2.74), and has the explicit form (2.82). The operator V{t) obeys 
the evolution equation (2.85), where the coupling Hamiltonian JC{t) is defined by Eq. (2.86). The 
matrix (2.87) generalizes according to Eq. (3.60): 



K = 


\afQ<^G<!>*Q aQ^G^Q* 
{afQ*^*G^*Q \a\^Q*^*G<^Q* 



|a| 


K{t) 


(Qio) 0 \fuUt) 0 \ 

VO Q*{0)) \ 0 U^{t)J 

^ I 0 \/q(0) 0 

"" \ia*f^*it)G<^*{t) \a\^<P*it)G<^{t)j V 0 1 0 Q*{0) 


(3.61) 


In terms of a complete set of single-particle states orthogonal to the condensate mode, the coupling 
Hamiltonian takes on the explicit form. 


}C{t) = \af a]ak{xj{t)\^{t)G<!>*{t)\xk{t)) 
j,k>l 

{i(^*)^^j^k{x*j{t)\^*{t)G^*{t)\xk{t)) + a^a]al{Xj{t)\^{t)G^{t)\xUt))'^ ■ 

j,k>l 

(3.62) 


B. Spin squeezing in the Bogoliubov approximation 

Despite the formal similarity of the single- and two-component cases, there is an important 
difference, which involves the special orthogonal mode of Eq. (3.2). In the case of most 

interest, when there is no single-particle coupling between the internal levels, i.e. (jJi 2 = 0 J 21 = 0, 
which we specialize to throughout this subsection, we can assume, as we discussed in Sec. HI A, that 
ai and a 2 are constants in time, and thus the mode |0) of Eq. (3.2) satisfies the GP equation (3.41). 
This allows us to make \<p{t)) one of the time-dependent modes orthogonal to the condensate 
mode Ixo(O) = it is convenient to choose |xi(t)) = If we further neglect the 

coupling of \(j)) to the other orthogonal modes, the coupling Hamiltonian (3.62) reduces to 

^{t) = \a\^fj{t)a\ai + 


(3.63) 

(3.64) 
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where we recall that = ai(0) = a^(o)) where we introduce the coupling parameter 
fj{t) = 

(3.65) 

= (m*{t)\m 

2 2 

= J (9ii\4>i\ +<?22l</’2l \(t>2\ 'jd^i (3.66) 

with 9 = arg(a!ia 2 /ct^) being a constant phase angle. The form (3.64), where we let ae*^ = |a|e*^, 
is the Hamiltonian of a free particle with momentum quadrature (aie~*^ + o|e*^)/\/2 and a variable 
mass; this Hamiltonian produces shearing and squeezing in the direction of the position quadrature 
at a variable rate given by 2 |q;| i/(t). 

We can solve for the Heisenberg-picture evolution of a^, 

V\t)a,V{t) = a,[l-i\afm] -ia\\afe^'>^m = . (3.67) 

In the first form, 

^{t) = \J^ 9{t') dt' (3.68) 

is the dimensionless integral of the coupling parameter. In the second form, 

S(C) = exp Q[C*a? - C(ai)^]) (3-69) 

is the single-mode squeeze operator [37-41], with the complex squeezing parameter given by ({t) = 
tj[t)e e ^ , where 

sinh 7 (t) = \af^{t) = tan[2z^(t)] . (3.70) 

These results determine the evolution operator V{t) up to a phase. Irrelevant though it is, the 
phase can be determined by integrating directly the evolution equation for V{t) or by considering 
the vacuum expectation value of V{t), with the result that 

V(t) = , (3.71) 

where v{t) = 5 |a|^C(^) — ^ tan~^[|a|^^(t)]. 

Suppose now that the orthogonal mode \(j)) begins in vacuum, so that the initial extended 
catalytic state is |'0ecs(6)) = !«) '?^(0))o® |vac, (^(0))i. Then the extended catalytic state at time t is 

l'0ecs(i)) =^gp(^)V(t)|i/’ecs(0)) = 6*"^*VX*))o ® ^gp±(O'?(C(i)) I vac, <^(0))i . (3.72) 

where 

oo 

^gp±(^) = . (3.73) 

n =0 

The physical state at time t follows from projecting onto the Al-particle sector, as specified by 
Eq. (2.8), giving 

l^7v(i)) = Y1 i(^>^(O)|5'(C(t))|vac,0(O))i (3.74) 
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The matrix elements of the squeeze operator can be evaluated explicitly (only even values of M have 
nonzero matrix elements), but we do not bother with that here, noting instead that the Bogliubov 
approximation requires that these matrix elements be small for M > 2. A more quantitative 
statement is that these matrix elements must be < l/\/iV for M > 2, which translates to small 
squeezing with | 7 (f)| < l/\/iV or In practice, since the elements of G are nearly 

equal (they typically differ only by small differences in scattering length for the hyperfine levels), 
the smallness of or \fi{t)\ is governed by the difference in the probability densities for the two 
internal levels, |(?!)i(x)|^ and |(/) 2 (x)|^; roughly speaking, the Bogliobov approximation requires that 
the two hyperfine levels not be separated spatially. 

If |C(f)| becomes too large, perhaps due to spatial separation of the internal levels, one can do 
a better job by returning to the Hamiltonian (3.63) and recalling that it arises, in the Bogoliubov 
approximation, from replacing and by a and a* in the original Schrodinger-picture Hamil¬ 
tonian. Restoring, in normal order, the creation and annihilation operators for the condensate 
mode to /C(t) gives a Kerr-like interaction between the condensate mode \(p) and the orthogonal 
mode \4>), 


m 


\n(t) 

\n(t) 


/oft I 210/ t\2 2 I -2ie, t^2 2\ 

( 2a),a|aoai-I-e (aj) Oq + e (a),) aij 

/ 10 t I —id t \2 !_/ s/ f I t \ 

(e a'oo + e a^ai) --??(t)(a(,ao + ajai) , 


(3.75) 

(3.76) 


where Qq = a^(o)- ^i^st term in Eq. (3.75) comes from scattering of \(j))- and |(^)-particles off 
one another, the second term from scattering of two |(/>)-particles into the |(^)-mode, and the last 
term from scattering of two |(^)-particles into the |(?!))-mode. 

The Hamiltonian iC(t) conserves the total particle number M = o^Qq + a\ai. Thus, if the initial 
state is in the A^-particle sector of the two modes, it stays there, and we can omit the projection 
onto the A^-particle sector that is involved in our use of an extended catalytic state. Moreover, the 
last term in Eq. (3.76) is proportional to M and thus becomes the c-number Nfi{t)/2] this term 
only introduces an irrelevant overall phase, so we can neglect it. Einally, we let a, and 02 be 
real and positive, which makes 0 = 0, thus leaving us with the Hamiltonian 


where 


^ss(i) = , 

•^x= 2 (®0®1 + “l®o) 


(3.77) 


(3.78) 


is the x-component of the Schwinger pseudo-spin of the two modes. The other two Schwinger 
operators are Jy = — f(aQai — ajag)/2 and = (agOg — OiOi)/2. 

At this point the model has been reduced to two modes, both of which participate in the dynam¬ 
ics. The single-particle states for the two modes, \4>{t)) and |^(t)), change in time according to the 
GP equations (3.35) and (3.41), both of which are expressions of the coupled GP equations (3.30) 
(with UJ 12 = ijJ 2 i = 0). The state of the two modes changes according to the evolution operator 
V{t) of Eq. (2.85), where one uses the two-mode Hamiltonian /Css(t), which has fixed creation and 
annihilation operators; this evolution is followed by application of the operator Z7gp(t) of Eq. (2.82), 
which translates the two-mode state to the modes that apply at time t. 

This model ignores the coupling of the two dominant modes to the other orthogonal modes, 
but its chief problem lies in an inconsistency between the GP equations (3.30) and the evolution 
under the Hamiltonian (3.77). The GP equations contain the quantities |ai| and |a 2 l > which can 
be interpreted as mean particle numbers for the two internal levels, Ni = la^l^ and N 2 = |a 2 pj 
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the internal levels initially have number uncertainties of order \/N^ and \/N 2 , which are small 
compared to the mean particle numbers. As discussed previously, the GP equations (3.30) leave 
\ai\^ and |a2|^ unchanged when lui 2 = = 0. The inconsistency arises because the number 

operators for the two internal levels are generally not conserved by For the model to be 

consistent, the number operators for internal levels, Afi = b\bi and A/’2 = ^2^2; should be conserved 
or nearly so, so that the mean particle numbers don’t change and the uncertainties remain small. 

Letting cti/a = cos{^p/2) and 02/0 = sin((/9/2), we have, from Eqs. (3.10) and (3.11), bi = 
cos(<y9/2) + ai sin((/j/2) and 62 = 62,02 ~ sm{ip/2) — cos{ip/2) and thus 

Ni = 6|6i = cos (^ + Ja; sin (^ , 

^ (3.79) 

M 2 = 6362 = -AA — cos Lp — sin (p . 

Thus the condition for the model to be consistent is that the condensate mode be an equal su¬ 
perposition of the two internal levels, i.e., p = 7r/2, or nearly so. The combination of the GP 
equations (3.30) for evolving the spatial mode functions and the Hamiltonian (3.77) to evolve 
the two-mode state in the case of an equal superposition of the internal levels is called the two- 
component formalism (or two-mode approximation) [12, 42, 43]. The two-mode approximation is 
more robust than the Bogoliubov-approximation squeezing results summarized in Eq. (3.74). In 
the context of condensates isolated in fairly well separated trapping potentials, very recent work 
has analyzed the effect of including four macroscopically excited modes, two in each well [44]. 

The term in is the so-called one-axis-twisting Hamiltonian [45]; it induces spin squeezing 
in states that are initially maximally polarized along the spin z axis, as is the case for an initial 
state that has all N particles in the condensate mode \^). The one-axis-twisting Hamiltonian is 
widely used to generate spin squeezing in BEGs [33-36] . Of particular relevance to our formulation 
are analyses of the interplay of spatial and spin dynamics [46-48]; in addition, Sinatra et al. [49] 
showed that the amount of squeezing is bounded from above by the initial noncondensed fraction 
at finite temperature. 


IV. SUMMARY AND CONCLUSION 

In this paper we develop a new framework for deriving the number-conserving Bogoliubov ap¬ 
proximation for a dilute-gas BEG. Our approach begins by introducing the extended catalytic 
state (2.4), a coherent state for the condensate mode and an arbitrary state for the modes or¬ 
thogonal to the condensate mode. The physical state with exactly N particles is retrieved from 
the extended catalytic state by projecting into the V-particle sector, as in Eq. (2.5). To formu¬ 
late the Bogoliubov approximation, we introduce the time-dependent interaction picture (2.18) 
in which the condensate mode is displaced to the vacuum. The field operators are thus of order 
and we can organize the BEG Hamiltonian in powers of . Requiring the terms of order 

to vanish yields the Gross-Pitaevskii equation (2.33). Going to the next order, , gives the 
conventional Bogoliubov Hamiltonian of Eqs. (2.36) and (2.37). Introducing the auxiliary Hamilto¬ 
nian (2.39) removes the unwanted phase diffusion from the conventional Bogoliubov Hamiltonian 
without affecting the physical state in the V-particle sector (to order N^). The result is the 
number-conserving Bogoliubov Hamiltonian of Eqs. (2.54) and (2.55). Analysis of the dynamics 
under the number-conserving Bogoliubov Hamiltonian shows that its GP part (2.56) evolves the 
single-particle, spatial mode structure forward in time, and the remaining part (2.58) evolves the 
state of these modes. The result is the particularly simple form (2.91) for the Schrodinger-picture 
evolution at Bogoliubov order. 
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In Sec. Ill we extend our approach to BECs with two internal levels. This turns out to be 
largely a matter of using a spinor notation that puts the derivation into a form that mimics the 
single-component derivation of Sec. II. Because of this formal similarity, generalization to multiple 
hyperfine levels would be straightforward. In Sec. Ill B we specialize the Bogoliubov approximation 
to the two dominant modes, the condensate mode and the mode orthogonal to it in the same two- 
dimensional subspace of the internal levels. We discuss how to generalize beyond the Bogoliubov 
approximation to the two-mode approximation for these two dominant modes, thus allowing a 
treatment of the spin squeezing of these two modes. 

The number-conserving Bogoliubov Hamiltonian we find is identical to the one derived originally 
by Castin and Dum [10]. The chief difference between our approach and that of Castin and 
Dum is that they worked in the Heisenberg picture, whereas our derivation is carried out in the 
Schrodinger picture and a closely allied interaction picture. There are several reasons for presenting 
a new framework for a derivation of the same result. The first is the modest one that the new 
derivation might highlight assumptions from a different perspective and ease the way forward on 
different, but related problems. Indeed, in our approach, the Bogoliubov excitations are easily seen 
as excitations on top of a pure condensate mode that evolves according to the GP equation. The 
method we nse for handling the phase diffusion that arises from assigning a phase to the condensate 
mode comes directly from re-asserting the number conservation that applies to a lossless BEG. A 
second reason is that working in the Schrodinger picture allows us to separate cleanly, within 
the Bogoliubov approximation, the evolution of the spatial mode structure from the evolution 
of the state of these modes. A central problem of dilute-gas BEG theory and of many other 
problems in many-body physics is how to do this separation appropriately, and our analysis can 
be instrnctive in how to formulate this separation. A third reason is that the Schrodinger picture 
allows us to identify the entanglement between the condensate mode and the orthogonal modes 
that are excited by the Bogolinbov Hamiltonian. In our approach, this entanglement arises when 
the extended catalytic state is projected into the Wparticle sector to obtain the physical state of 
the BEG. The entanglement is something that can be analyzed easily in our Schrodinger-picture 
formulation, whereas even the proper formulation of entanglement is difficnlt in the Heisenberg 
picture. 

These three reasons motivated our work on this topic. We trust that they justify its presentation 
to the wider scientific community. 


ACKNOWLEDGMENTS 

The anthors thanks A. B. Tacla for useful and stimulating conversations. This work was sup¬ 
ported in part by National Science Eoundation Grants No. PHY-1212445 and No. PHY-1314763 
and by Office of Naval Research Grant No. N00014-15-1-2167. 


[1] N. N. Bogoliubov, On the theory of superfluidity, Journal of Physics (USSR) 11, 23 (1947). 

[2] A. L. Fetter, Nonuniform states of an imperfect Bose gas, Annals of Physics 70, 67 (1972). 

[3] K. Huang, Statistical Mechanics, 2nd Edition (Wiley, 1987). 

[4] Y. Nambu, Quasi-particles and gauge invariance in the theory of superconductivity. Physical Review 
117, 648 (1960). 

[5] J. Goldstone, Field theories with superconductor solutions, II Nuovo Cimento 19, 154 (1961). 

[6] M. Girardeau and R. Arnowitt, Theory of many-boson systems: Pair theory, Physical Review 113, 755 
(1959). 



28 


[7] C. W. Gardiner, Particle-number-conserving Bogoliubov method which demonstrates the validity of the 
time-dependent Gross-Pitaevskii eguation for a highly condensed Bose gas. Physical Review A 56, 1414 
(1997). 

[8] C. W. Gardiner, P. Zoller, R. J. Ballagh, and M. J. Davis, Kinetics of Bose-Einstein condensation in 
a trap, Physical Review Letters 79, 1793 (1997). 

[9] Y. Gastin and R. Dum, Instability and depletion of an excited Bose-Einstein condensate in a trap. 
Physical Review Letters 79, 3553 (1997). 

[10] Y. Gastin and R. Dum, Low-temperature Bose-Einstein condensates in time-dependent traps: Beyond 
the U(1) symmetry-breaking approach, Physical Review A 57, 3008 (1998). 

[11] Y. Gastin, Bose-Einstein condensates in atomic gases: Simple theoretical results, in Coherent Atomic 
Matter Waves, Les Houches - Ecole d’Ete de Physique Theorique No. 72, edited by R. Kaiser, G. West¬ 
brook, and F. David (Springer, Berlin, 2001) pp. 1-136. 

[12] A. S. Sprensen, Bogoliubov theory of entanglement in a Bose-Einstein condensate. Physical Review A 
65, 043610 (2002). 

[13] S. A. Gardiner and S. A. Morgan, Number-conserving approach to a minimal self-consistent treatment 
of condensate and noncondensate dynamics in a degenerate Bose gas. Physical Review A 75, 043621 
(2007). 

[14] S. A. Gardiner and T. P. Billam, Number-conserving approaches for atomic Bose-Einstein condensates: 
An overview, in Quantum Gases, Gold Atoms, Vol. 1 (Imperial College Press, 2011) pp. 133-145. 

[15] P. Mason and S. A. Gardiner, Number-conserving approaches to n-component Bose-Einstein conden¬ 
sates, Physical Review A 89, 043617 (2014). 

[16] M. J. Steel, M. K. Olsen, L. 1. Plimak, P. D. Drummond, S. M. Tan, M. J. Collett, D. F. Walls, and 
R. Graham, Dynamical quantum noise in trapped Bose-Einstein condensates. Physical Review A 58, 
4824 (1998). 

[17] A. Sinatra, C. Lobo, and Y. Gastin, The truncated Wigner method for Bose-condensed gases: limits 
of validity and applications. Journal of Physics B: Atomic, Molecular and Optical Physics 35, 3599 
( 2002 ). 

[18] A. J. Leggett, Bose-Einstein condensation in the alkali gases: Some fundamental concepts. Reviews of 
Modern Physics 73, 307 (2001). 

[19] A. J. Leggett, The relation between the Gross-Pitaevskii and Bogoliubov descriptions of a dilute Bose 
gas. New Journal of Physics 5, 103 (2003). 

[20] J. Dziarmaga and K. Sacha, Bogoliubov theory of a Bose-Einstein condensate in the particle represen¬ 
tation, Physical Review A 67, 033608 (2003). 

[21] J. Dziarmaga and K. Sacha, Images of a Bose-Einstein condensate: diagonal dynamical Bogoliubov 
vacuum. Journal of Physics B: Atomic, Molecular and Optical Physics 39, 57 (2006). 

[22] Z. Jiang, A. B. Tacla, and C. M. Caves, Bosonic particle-correlated states: A nonperturbative treatment 
beyond mean field, in preparation. 

[23] P. Villain, M. Lewenstein, R. Dum, Y. Gastin, L. You, A. Imamoglu, and T. A. B. Kennedy, Quantum 
dynamics of the phase of a Bose-Einstein condensate, Journal of Modern Optics 44, 1775 (1997). 

[24] 1. Danshita, K. Egawa, N. Yokoshi, and S. Kurihara, Gollective excitations of Bose-Einstein condensates 
in a double-well potential. Journal of the Physical Society of Japan 74, 3179 (2005). 

[25] F. Trimborn, D. Witthaut, and H. J. Korsch, Exact number-conserving phase-space dynamics of the 
M-site Bose-Hubbard model. Physical Review A 77, 043631 (2008). 

[26] F. Trimborn, D. Witthaut, and H. J. Korsch, Beyond mean-field dynamics of small Bose-Hubbard 
systems based on the number-conserving phase-space approach. Physical Review A 79, 013608 (2009). 

[27] B. Oles and K. Sacha, N-conserving Bogoliubov vacuum of a two-component Bose-Einstein condensate: 
Density fluctuations close to a phase-separation condition. Journal of Physics A: Mathematical and 
Theoretical 41, 145005 (2008). 

[28] J. Schachenmayer, A. J. Daley, and P. Zoller, Atomic matter-wave revivals with definite atom number 
in an optical lattice. Physical Review A 83, 043614 (2011). 

[29] T. P. Billam, P. Mason, and S. A. Gardiner, Second-order number-conserving description of nonequilib¬ 
rium dynamics in finite-temperature Bose-Einstein condensates. Physical Review A 87, 033628 (2013). 

[30] Z. Jiang, Particle Gorrelations in Bose-Einstein Gondensates, Ph.D. thesis. University of New Mexico 
(2014). 

[31] L. Pitaevskii and S. Stringari, Bose-Einstein Gondensation (Oxford University Press, Oxford, 2003). 


29 


[321 M. Lewenstein and L. You, Quantum phase diffusion of a Bose-Einstein condensate. Physical Review 
Letters 77, 3489 (1996). 

[33] J. Esteve, C. Gross, A. Weller, S. Giovanazzi, and M. K. Oberthaler, Squeezing and entanglement in 
a Bose-Einstein condensate, Nature 455, 1216 (2008). 

[34] J. Grond, J. Schmiedmayer, and U. Hohenester, Optimizing number squeezing when splitting a meso¬ 
scopic condensate, Physical Review A 79, 021603 (2009). 

[35] M. F. Riedel, P. Bohi, Y. Li, T. W. Hansch, A. Sinatra, and P. Treutlein, Atom-chip-based generation 
of entanglement for quantum metrology. Nature 464, 1170 (2010). 

[36] C. Gross, T. Zibold, E. Nicklas, J. Esteve, and M. K. Oberthaler, Nonlinear atom interferometer 
surpasses classical precision limit. Nature 464, 1165 (2010). 

[37] A. M. Perelomov, Generalized coherent states and some of their applications, Soviet Physics Uspekhi 
20, 703 (1977). 

[38] J. N. Hollenhorst, Quantum limits on resonant-mass gravitational-radiation detectors. Physical Review 
D 19, 1669 (1979). 

[39] C. M. Gaves and B. L. Schumaker, New formalism for two-photon quantum optics. I. Quadrature 
phases and squeezed states. Physical Review A 31, 3068 (1985). 

[40] B. L. Schumaker and G. M. Gaves, New formalism for two-photon quantum optics. II. Mathematical 
foundation and compact notation. Physical Review A 31, 3093 (1985). 

[41] B. L. Schumaker, Quantum mechanical pure states with Gaussian wave functions. Physics Reports 135, 
317 (1986). 

[42] U. V. Poulsen and K. Mplmer, Positive-p simulations of spin squeezing in a two-component Bose 
condensate, Physical Review A 64, 013616 (2001). 

[43] U. V. Poulsen, Bose-Einstein Gondensates: Excursions Beyond the Mean Field, Ph.D. thesis. University 
of Aarhus (2002). 

[44] J. Gillet, M. A. Garcia-March, T. Busch, and F. Sols, Tunneling, self-trapping , and manipulation of 
higher modes of a Bose-Einstein condensate in a double well. Physical Review A 89, 023614 (2014). 

[45] M. Kitagawa and M. Ueda, Squeezed spin states. Physical Review A 47, 5138 (1993). 

[46] A. Sinatra and Y. Gastin, Binary mixtures of Bose-Einstein condensates: Phase dynamics and spatial 
dynamics. The European Physical Journal D 8, 319 (2000). 

[47] A. Sprensen, L.-M. Duan, J. 1. Cirac, and P. Zoller, Many-particle entanglement with Bose-Einstein 
condensates, Nature 409, 63 (2001). 

[48] Y. Li, P. Treutlein, J. Reichel, and A. Sinatra, Spin squeezing in a bimodal condensate: spatial dynamics 
and particle losses. The European Physical Journal B 68, 365 (2009). 

[49] A. Sinatra, E. Witkowska, J.-G. Dornstetter, Y. Li, and Y. Gastin, Limit of spin squeezing in finite- 
temperature Bose-Einstein condensates. Physical Review Letters 107, 060404 (2011). 


